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Density functional theory calculations of NdAl0 3 /SrTi0 3 heterostructure show that two-dimensional 
electron gas (2-DEG) is produced at the interface with a built-in potential of ~0.3 eV per unit cell. The 
effects of surface defects on the phase stability and electric field of 2-DEG have been investigated. It is found 
that oxygen vacancy is easily to form on the NdAlO 3 (001) surface, with a low threshold displacement energy 
and a low formation energy. This point defect results in surface reconstruction and the formation of a zigzag 
-Al-O-Al- chain, which quenches the built-in potential and enhances the carrier density significantly. These 
results will provide fundamental insights into understanding how surface defects influence the electronic 
behavior of 2-DEG and tuning their electronic properties through surface modification. 

With the breakthrough of cutting-edge thin-film epitaxial techniques like pulsed laser deposition and 
molecular-beam epitaxy, high quality heterointerfaces based on perovskite transition metal oxides 
have been achieved experimentally and revealed to exhibit a large number of fascinating electrical and 
magnetic properties 1 . A prominent example is the discovery of the two-dimensional electron gas (2-DEG) at the 
interface between two band insulators, polar LaAl0 3 (LAO) and non-polar SrTi0 3 (STO) 2 . Due to its intriguing 
functionalities 3 " 5 , unparalleled environmental stabilities over traditional semiconductors, and capability for 
nanoscale read/write/erase 6 , as well as field-effect transistor and non-volatile memory effects, oxide heterostruc- 
tures are now considered as promising candidates for breeding next- generation high-performance electronics 
and spintronics. 

In recent years, the origin of the 2-DEG formation at the LAO/STO(001) interface has attracted broad 
interests 7 ' 8 . Along the [001] direction STO is composed of neutral layers (SrO)° and (Ti0 2 ) 0 , and LAO consists 
of alternately charged atomic layers (LaO) + and (A10 2 )~. Electronic reconstruction can occur near the electron- 
rich n-type LaO/Ti0 2 interface to avoid the polarization catastrophe caused by the diverging Coulomb field 9-12 . 
On the other hand, oxygen vacancies created during the deposition process can dope electrons in STO substrate 
and induce the interface conductivity 1314 . 

Cation nonstoichiometry and intermixing at the interface 15 " 18 have also been proposed to result in the con- 
ducting interface. To investigate the role of the AB0 3 polar layer and explore 2-DEG with higher carrier density 
and hole mobility, various perovskite oxide thin films on STO substrate, e.g., NdAl0 3 19 , GdTi0 3 20 , LaMn0 3 21 , 
LaGa0 3 22 ' 23 , and LaCr0 3 24 , have been studied. Recently, Annadi et al reported that NdAl0 3 /SrTi0 3 (NAO/STO) 
can also produce 2-DEG at the interface that exhibits a stronger localization of 2-DEG than LAO/STO inter- 
faces 19,25 . However, there is a lack of theoretical understanding of the formation mechanism. 

In the potential application of 2-DEG as electronic devices in aerospace industries, surface defects may be 
created due to radiation in a harsh environment. Knowledge of radiation response, surface defect stability and its 
influence on 2-DEG behavior are thus of vital importance for realistic applications. This is also important for 
surface modification by ion beam implantation 26 to improve overall performance of future electronic devices 
under harsh environments and to design novel materials with enhanced performance. In spite of the studies on 
the stabilities of the intrinsic defects on surface or at interfaces 27 " 30 , very few experimental and theoretical 
investigations of heterostructures under irradiation have been reported. In this work, a first-principles study 
of the NdAlO 3 /SrTiO 3 (001) heterointerface is carried out to explore its low energy radiation response and the 
stability of the built-in potential with respect to surface defect formation. Our results reveal that the surface 
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oxygen within the NAO/STO heterostructure can be easily displaced, 
which reduces the built-in potential to zero, and enhances the carrier 
density significantly. However, the surface aluminum vacancy 
formation is hard to form, and its formation energy is relatively high, 
independent of the thickness of NAO, and is shown to sustain the 
built-in potential. 

The first-principles calculations are performed within the density 
functional theory (DFT) framework using the projector augmented 
wave method, as implemented in the Vienna Ab Inito Simulation 
Package (VASP) 31 . The exchange- correlation effects are treated by 
the generalized- gradient approximation (GGA) in the Perdew- 
Burke-Ernzerhof parametrization 32 , with spin-polarized effects con- 
sidered. A plane-wave cutoff energy of 400 eV and a 4 X 4 X 1 
k-point sampling in Brillouin zone are used. Slab models for 
(NAO) n /(STO) 4 (n=2, 4, 6, 8, and 10), in which one slab consisting 
of n layers of NAO and four layers of STO, and a vacuum region of 
15 A separating the periodically repeated slabs, are employed in the 
calculations. A y/l x \fl unit cell is taken into account. We only 
consider the n-type interface, which is formed by a layer sequence 
of . . .-(SrO)-(Ti0 2 )-(NdO)-(A10 2 )-. . .. The surface of the NAO over- 
layer is A10 2 terminated. The geometry of the simulation cell for 
(NAO) 4 /(STO) 4 is shown in Fig. 1(a). The two-dimensional lattice 
parameters are fixed at 3.94 A, which is the optimized lattice con- 
stant for the bulk STO. Structural optimization is carried out in the 
condition that the Hellmann-Feynman force on each atom is smaller 
than 0.01 eV A" 1 . The dipole correction is applied to remove the 
undesired electrostatic interaction between periodic cells along the 
direction perpendicular to the surface. Correlation effect 33 has been 
tested and it is shown that the obtained results are very similar to 
those without considering this effect. 

The total density of states for 2-4 layers of NAO film on STO(001) 
presented in Fig. 1(b) show that (NAO) 2 /(STO) 4 owns a band gap as 
small as 0.6 eV. The band gap, defined by occupied O 2p states at the 
surface and unoccupied Ti 3d states at the interface, decreases with 
increasing thickness of the NAO film and a crossover to an insulator- 




Figure 1 | (a) Structural model for (NAO) 4 /(STO) 4 and (b) total density 
of states of (NAO) n /(STO) 4 , where the Fermi level E F is set to zero. 



to -metal transition occurs for thin film thickness between 3 and 4 
unit cells (u.c). Experimentally it is also observed that a transition 
from insulating to metallic interface occurs at a thickness around 4 
u.c 19 . Such agreement is surprising since our GGA calculations 
underestimate the band gap of STO by 1.95 eV and the critical thick- 
ness should increase by around two layers when the valence-band 
offset at the interface is correctly described within GGA and the true 
energy gap of STO is considered 34 . There is a possibility that the true 
valence-band offset is very different from the GGA results such that 
the energy gap underestimation is effectively cancelled by the error in 
the valence-band offset 34 . Layer-resolved density of states for 
(NAO) 4 /(STO) 4 interface is presented in Fig. 2(a). It is shown that 
the valence edge in NAO, which is mainly contributed by O 2p in the 
surface A10 2 layer, increases linearly in energy with the distance from 
the interface. This built-in potential (—0.3 eV/u.c.) results in the 
emergence of a small concentration of charge carrier, i.e., holes at 
the top NAO layer and electrons that spread over several STO layers, 
supporting the idea that electronic reconstruction could be the driv- 
ing mechanism for the creation of the 2-DEG at the interfaces 9 " 11 . It is 
found that NAO/STO heterointerface exhibits similar characteristics 
in the insulator-metal transition and the thickness dependency of the 
polar layer to LAO/STO system 10,34 . 

To investigate the low energy radiation response of NAO/STO 
heterostructure, the threshold displacement energies (E d ) of surface 
atoms are calculated. Threshold displacement energy is defined as 
the minimum amount of transferred kinetic energy necessary to 
permanently displace an atom from its original lattice site, thus, 
forming a stable defect 35 . This physical parameter is important for 
determining the primary state of radiation damage in a material and 
estimating the total amount of damage states created under ion irra- 
diation 36 " 38 . Along the direction perpendicular to the surface, the E d 
are calculated to be 9 eV for surface oxygen and 40 eV for surface 
aluminum, which correspond to electron irradiation of 60 and 
360 KeV, respectively. The respective created defects are mainly sur- 
face oxygen vacancy and aluminium vacancy. The low E d value of 
9 eV for surface oxygen indicates that it is easily to be displaced and 
surface oxygen vacancies will be one of the dominant defects under 
irradiation. In order to study the stability of surface defects in NAO/ 
STO heterostructure, one oxygen or aluminum is removed from the 
surface, and the formation energies for surface oxygen vacancy and 
aluminum vacancy are calculated. The formation energy is defined 
by Ef = E de f — E per f + /Hi 29 . Here E de f and E per f are the total energies for 
defective and perfect heterostructures, respectively. The \i { is the 
chemical potential of oxygen or aluminum, which is obtained under 
O-rich condition 39 . The defect formation energies as a function of 
NAO layers are provided in Fig. 3(a). It should be pointed out that the 
variation of defect formation energy with thin film thickness is inde- 
pendent of the oxygen or Al chemical potential. As shown in the 
figure, the formation of surface oxygen vacancy becomes easier with 
increasing NAO thin film thickness. This is consistent with the 
theoretical calculations of LAO/STO performed by Zhang et al. 29 
However, the formation of surface Al vacancy shows negligible 
dependence on the thin film thickness. It is also found that the 
formation energy for surface Al vacancy is much higher than that 
for surface O vacancy, confirming the finding that surface oxygen 
vacancies will be one of the dominant defects under irradiation. 

The vertical displacement in the NAO layers for an ideal (NAO) 4 / 
(STO) 4 and the heterointerfaces with surface oxygen or aluminum 
defects is investigated in Fig. 3(b) and 3(c). For the ideal (NAO) 4 / 
(STO) 4 structure, it is noted that Nd relaxes outward by 0.10-0.22 A 
and the Nd far from the interface shows the largest relaxation of 
0.22 A. On the other hand, the displacement of Al is small and the 
surface Al only shifts inward by 0.05 A. Oxygen ions in NdO and 
A10 2 layers relaxes inward by 0.01-0.05 A and 0.02-0.12 A, respect- 
ively. This large polar distortion has a remarkable impact on the 
electronic behavior of the heterostructures and results in electronic 
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Figure 2 | Layer- resolved density of states for (a) ideal (NAO) 4 /(STO) 4 interface; (b) interface with surface oxygen vacancy; (c) interface with surface 
aluminum vacancy, where the Fermi level E F is set to zero. 



reconstruction. The polar-type distortion in NdO layers is very sim- 
ilar to that in LaO layers of LAO/STO heterointerfaces 10 . The intro- 
duction of surface O vacancy is shown to have slight influences on the 
vertical displacement of ions in NdO and A10 2 layers. When oxygen 
vacancy exists on the surface, NdO layers show small relaxation of 



0.04-0.06 A. The surface and subsurface A10 2 layers are found to 
exhibit different features, i.e., surface A10 2 layer relaxes inward by 
0.07 A, while subsurface A10 2 layers relax outward by 0.05 A. The 
most striking result for NAO/STO heterostructure with surface O 
vacancy is that oxygen ions in NdO layers have large in-plane dis- 
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Figure 3 | (a) Variation of vacancy formation energy with thin film thickness; (b) vertical displacements of cations in the NAO layers for (NAO) 4 / 
(STO) 4 , and (c) vertical displacements of anions in the NAO layers for (NAO) 4 /(STO) 4 . 
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Figure 4 | Illustration of optimized geometrical configuration for 
NdA10 3 thin film with surface oxygen vacancy, where the atom 
representations are the same as those in Fig. 1(a). 

placements. According to Fig. 4, the oxygen ions in the NdO layers 
are displaced by 0.40 A from their ideal lattice sites. This displace- 
ment exhibits a compression and expansion behavior alternately, 
forming a zigzag -Al-O-Al- chain structure which has also been 
observed experimentally in LAO/STO interface 40 . Cen et al. have 
performed first-principles density functional theory calculations of 
LAO films on STO substrates, and a zigzag -Al-O-Al chian structure 
is also found when oxygen ions are removed from the LAO surface 
layer 41 . In addition, the large distortion is expected to significantly 
affect the built-in potential and electronic properties of heterostruc- 
tures. For a NAO/STO heterostructure with surface Al vacancy, the 
optimized structure shows nearly zero in-plane displacement. This 
phenomenon is different from the case of surface O vacancy, but it is 
very similar to that of an ideal NAO/STO. Vertically, the ions in 
different NAO overlayers are also found to relax in a similar way 
to that in ideal NAO/STO. Nevertheless, the Nd ions farthest from 
the interface are found to have larger outward relaxation of 0.59 A, as 
compared with the value of 0.33 A in an ideal NAO/STO structure. 
Another difference is that the surface Al shifts outward and inward in 
the defective and ideal heterostructure, respectively. The electronic 




Number of NAO layers 

Figure 5 | Variation of carrier densities with thin film thickness for ideal 
and defective NAO/STO. The O vac and Al vac represent oxygen and 
aluminum vacancy, respectively. 



properties for NAO/STO with surface Al vacancy are thus expected 
to be similar to those for ideal states. 

To explore the built-in potential stability with respect to vacancy 
defects, the layer-resolved density of states for defective heterostruc- 
ture are compared with those for ideal states in Fig. 2. It is of interest 
to find that the surface oxygen vacancy causes significant electron 
redistribution and the built-in potential vanishes. This is because the 
formation of surface oxygen vacancy in the A10 2 layer introduces 
two excess electrons, which interact with the extra holes near the 
valence band maximum. The Coulomb coupling between the oxygen 
vacancy donor and the host results in charge transfer and hole com- 
pensation. This effect reduces the polarization of the NAO overlayer, 
indicating that the carrier density will be increased and the conduc- 
tivity at the interface will be enhanced. In contrast, a neutral Al 
vacancy at the A10 2 surface would produce three holes and leaving 
the overall charge unchanged. The tilt of the NAO band structure due 
to the built-in potential is thus sustained because the A10 2 surface 
charge cannot be compensated by Al vacancy. 

By adding the occupation numbers of the states corresponding to 
the electrons transferred from NAO to STO, the carrier density in 
STO can be obtained. Fig. 5 presents the variation of carrier densities 
with thin film thickness for ideal and defective NAO/STO. In all 
cases, it is found that the carrier density increases with the thin film 
thickness. Annadi et al. also observed that at low temperatures below 
20 K down to 2 K the carrier density increases as the NAO thickness 
increases from 6 to 12 u.c 25 . This is because with increasing NAO 
thickness the local valence band maximum of NAO in the surface 
region rises more above the conduction band minimum of STO, 
which results in larger charge transfer from NAO to STO and smaller 
electric filed inside NAO. For a 6 u.c. thin film, the carrier density is 
predicted to be 1.85 X 10 14 electrons/cm 2 , which is larger than the 
experimental value of 0.53 X 10 14 electrons/cm 2 . Similar phenom- 
enon has also been observed in LAO/STO system 34 . In the case of 
NAO/STO with surface oxygen vacancy, the carrier density is 
enhanced significantly as compared with the ideal states, due to 
the reduced built-in potential. Surface aluminum vacancy, on the 
other hand, lowers the density slightly because the introduced holes 
in the surface area lead to less electrons transferring from NAO to 
STO. These results suggest that the electronic behavior is influenced 
significantly under radiation environment and the surface of hetero- 
structure can be modified by introduction of defects. 

In summary, first-principles calculations show that two-dimen- 
sional electron gas can be produced at NAO/STO(001) heterointer- 
face with a critical thickness of 3 ~ 4 u.c. The threshold displacement 
energies for surface O and Al are determined to be 9 and 40 eV, 
respectively, and the formation energy for a surface Al vacancy is 
much higher than that for a surface O vacancy, suggesting that the 
heterointerface is susceptible to irradiation and oxygen vacancy is 
one of the dominant defects. Surface Al vacancy is found to have 
minimal effects on the built-in potential in NAO/STO, while surface 
oxygen vacancy reduces the potential to zero due to surface recon- 
struction and the formation of a zigzag -Al-O-Al- chain. The holes 
introduced by surface aluminum vacancy lead to less electrons trans- 
ferring from NAO to STO and lower the carrier density slightly; 
however, the excess electrons introduced by surface oxygen vacancy 
compensate the holes at the surface and enhance the carrier density 
significantly. These results provide an important understanding of 
radiation effects on the heterostructures, which is a necessity for 
achieving radiation-resistant electronic devices based on oxide 2- 
DEG through manipulating the electronic surface properties. 
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